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ABSTRACT 



A digital computer program for the geometrically 
nonlinear analysis of totally arbitrarily loaded 
shells of revolution (SATANS-II) was modified to more 
accurately account for the conditions at the pole of 
the shell. This program was used to determine the 
buckling load of shallow spherical shells of various 
sizes when subjected to static axisymmetric, dynamic 
axisymmetric, and dynamic nearly axisymmetric 
step- pressure loads of infinite duration. A 
comparison was made between the new buckling results 
and previous results obtained without the new pole 
routine. The comparison- revealed a significant change 
in the buckling pressures, due solely to the change in 
the pole routine. The new static axisymmetric, 
dynamic axisymmetric, and even the dynamic asymmetric 
critical buckling pressure loads appear to be fairly 
reliable results for perfect, shallow shells. 
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= nondimensional inplane stiffness 
= the modulus of elasticity of the shell 
= the rise of the spherical cap at the pole 
= the thickness of the shell 
= the mass density of the shell 

= the meridional bending moment per unit length 

= the Fourier index 

= a nondimensional applied load 

= the nondimensional critical pressure 

= the classical buckling pressure of a complete 

sphere 

= a column matrix containing the coefficients 
th 

of the n term in the series expansion of the 
applied load 

= the normal distance from the axis of revolution 
to the surface of the cap 
= the normal distance from the axis to the cap 

in the base plane; the maximum value of r 
= the radii of curvature in the s and0 

directions, respectively 
= the meridional distance along the surface 
of the shell 

= the nondimensional time 
- the time 
= a reference time 
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= the displacements in the s, Q and J 
directions, respectively 
= nondi me nsional series coefficients of U, V, H 
= a nondimensional measure of the volume of 
the shell deformation 
= the peak in the time history of the 

parameter V 

= the displacement in the S direction 
th 

in the n harmonic 
= the nondimensional time increment 
= distance between stations 

= the nondimensional parameter governing the 

magnitude of the load applied in the asymmetric 
harmonics 

= the coordinate normal to the surface of the shell 
= the circumferential angle measured about 
the axis of revolution 
= a nondimensional geometric parameter used 
to describe the spherical cap 
= Pcisson's ratio 

= the normal distance from the base plane to 
the middle surface of the undeformed cap 
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I. 



INTRODD CTION 



In 1973 a digital computer study was presented by Ball 
and Burt [1] for the dynamic buckling load of clamped 
shallow spherical shells subjected to axisymmetric and 
nearly axisymmetric step- pressure loads. A static buckling 
analysis of the same spherical shells had been carried out 
in 1970 by Stilwell and Ball [2]. In these two studies the 
digital computer program SATANS-I [3] was used to calculate 
the critical buckling pressures for a large range of shell 
sizes. Other studies of the buckling of shallow shells have 
been conducted by Huang [4,5], by Stephens and Fulton [6], 
by Lock et al. [7], by Stricklin [8], and most recently by 
Akkas [9]. In Reference 1 the results from these other 
studies, except for those by Akkas, are compared with the 
results from SATANS-I for both static and dynamic buckling. 
In the axisymmetric static analysis the comparison with the 
results obtained by Huang [4] revealed that the SATANS-I 
results were higher than Huang's results for several shell 
sizes. In the dynamic, axisymmetric buckling analysis the 
SATANS-I results again either agreed closely with, or were 
somewhat higher than , the results by Huang [5], Stephens 
and Fulton [6], and Stricklin [8]. However, it was noted 
then that there was a general lack of consistent agreement 
among any of the sets of results. As a consequence, it 
appeared at that time that the axisymmetric buckling problem 
had not yet been totally resolved and that additional 
studies would be appropriate. 

In the asymmetric dynamic buckling analysis of Reference 
1 the few comparisons that could be made for the critical 
load also indicated that the SATANS-I results may be too 
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high. A ccmparison of the recent estimates for the 
asymmetric dynamic buckling load obtained by Akkas [9] with 
the SATAMS-I results also reveals the SATANS-I results to be 
well above those of Akkas [9]. However, it should be noted 
that the results obtained by Akkas were from his attempt to 
obtain a lower bound on the critical asymmetric load. This 
bound on the buckling load is obtained without the execution 
of a complete transient response analysis on the asymmetric 
part of the response of the shell, as is done in SATANS-I. 
In Akkas* analysis (Problem 1) the transient nonlinear 
axisymmetric response is computed, and a determinant is 
examined for possible bifurcation into asymmetric motion at 
each time step. The minimum load at which the determinant 
becomes zero is defined as the lower bound of the critical 
load. 

As a consequence of the generally high buckling loads 
predicted by SATANS-I, a re-examination of the static and 
dynamic buckling of the shallow spherical shell was made in 
an attempt to determine the possible cause, or causes, of 
the high buckling loads. In our search we discovered that a 
modification of the manner in which the pole conditions are 
numerically approximated significantly lowered the buckling 
loads to values that are now in good agreement with the 
other results. The new procedure for handling the pole 
condition is given in section III of this thesis. The new 
buckling results are given in section V. 

In addition to the pole condition modifications and the 
new buckling results the author has also made another 
significant change to the SATANS family of codes. In 
particular, the SATANS-II program for the geometrically 
nonlinear analysis of totally arbitrarily loaded shells of 
revolution, developed by Hyan [10] in 1972 to handle more 
complex and larger problems, was modified to make the 
computer memory requirement a variable quantity. This 
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quantity is specified by the user to fit the particular 
problem being run. It eliminates the large core requirement 
of SATANS-II for small problems and allows for much larger 
problems to be solved than could be solved by SATANS-II. 
The new program with the pole condition and memory 
modifications will hereafter be called SATANS-IIA. It is 
described in section II. 
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II. ASCRIPTION OF SATAH S-IIA 



SATANS-II was developed by Ryan [10] from SATANS-I and 
incorporated the full trigonometric expansion of the applied 
load and solution vector, and introduced the handling of 
imperfections into the code. These modifications allow the 
analysis of shells under totally arbitrary loads, as well as 
imperfection studies on actual shells with measured 
imperfections [11]. Unfortunately, the original deck of 
cards for SATANS-II was destroyed. Professor Johann Arbocz 
of CALTECH had a listing of SATANS-II and punched a deck of 
cards with the changes to SATANS-I given in that listing. A 
copy of this deck was sent to Professor Ball. These cards 
have been added by the author to the original SATANS-I 
described by Ryan [10] and a complete version of SATANS-II 
has been reconstructed. SATANS-IIA is a modification by the 
author of the reconstructed SATANS-II program. A listing of 
SATANS-IIA can be found in Appendix A. The listing contains 
an example problem for the dynamic analysis of a clamped, 
truncated cone subjected to an impulsive loading which is 
uniform along the meridian and varies in a cosine 
distribution over one-half of the circumference. This 
problem is a sample problem suggested by the Lockeed 
Missiles and Space Corp. [12]. A condensed version of the 
output from the example problem is given in Appendix B. 
Input data preparation for SATANS-IIA can be found in 
Appendix C. The basic users manual, which includes 
preparation of input subroutines and the theory of the 
program, is contained in Reference 3, which can be obtained 
through COSMIC (M70-10098, LAR-10736) , or ASIAC [13]. A 
users manual which includes preparation and handling of 
imperfection data within the SATANS programs can be found in 
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Ref. [10]. The above information, along with the following 
discussion, will inform the user on the capabilities and 
proper use of SATANS-IIA. 

The modification of the SATANS-II program to make its 
core requirement variable was accomplished by putting in a 
single dimension statement at the beginning of the program, 
with subsequent dimensioning within the subroutines to only 
the first element of the vector or matrix. This is a 
convenient feature of the FORTRAN-IV language in which the 
program is written. The actual vector and matrix sizes are 
transmitted to the subroutines by ■ an individual parameter 
list. Construction of the initial dimension statement and 
core request size is as follows: 

The basic size of the program on the IBM-360/67 Digital 
Computer, without the initial dimension statement, is 
272,000 bytes. This figure includes approximately 19,000 
bytes of buffer space required for execution. Within the 
main dimension statement are fifteen variables. However, 
only three parameters are needed to specify the sizes of 
these fifteen variables. 

Let a= The number of stations along the meridian of the 
shell times the number of harmonics considered. 

Let b= a, plus two fictitious stations times the number 

of harmonics considered. 

Let c= The number of harmonics considered. 

The main dimension statement would then be constructed 

as , 

DIMENSION P(4,4,a), DEE(4,4,a), DST(4,4,a), X(4,a), 
PHIXB(a), PHITB (a) , Z(4,b), Z0(4,b), 

Z2(4,b), Z3(4,b), ZD0T(4,b), IS(99,c), 
JS(99,C), ID(99,c), JD(99,c) 
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The 99 's above limit the user to 99 harmonics in any one 
run and an unlimited number of meridional stations. The 
core requirement for the general case would be, 

272,000 + 216a + 80b ♦ 1584c = bytes of core required. 

For a sample calculation of the core requirements 
consider the example of a spherical cap with 40 stations 
along the meridian, and an asymmetric analysis With two 
harmonics. Therefore, 

a= 40 (stations) x 2 (harmonics) = 80 
b= 80 + 2 X 2 (harmonics) = 84 
c= 2 (harmonics) 

Thus, for the variables P, DBF, DST, 

3 X (4x4x80) = 3840 (words) x 4 = 15,360 bytes 

for the variable X, 

4 X (80) = 320 (words) x 4 = 1280 bytes 

for the varibles PHIXB, PHITB, 

2 X (80) = 160 (words) x 4 = 640 bytes 

for the varibles Z, ZO, Z2, Z3, ZDOT, 

5 X (4x84) = 1680 (words) x 4 = 6720 bytes 

lastly, for the varibles ID, JD, IS, JS, 

4 X (99x2) = 792 (words) x 4 = 3168 bytes 

Therefore, the total size of the main dimension 

statement would be 27,168 bytes. This figure would be 
rounded up to the nearest even thousand bytes, i^. 28,000 
bytes. Finally, the core requirement for this example 
problem would be 

272,000 + 28,000 = 300,000 bytes. 
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III. IMPBOVED POLE ROUTINE 



The SATANS code is based upon Sander's geometrically 

nonlinear equations under the conditions of small strains 

and moderately small rotations. The formulation is in four 

second order nonlinear partial differential equations in 

terms of U, V, H, and a , where 0, V, and H are the 

s 

meridional, circumferential and normal displacements 

respectively, and a is the meridional bending moment. The 

s 

nonlinear partial differential equations in the coordinates 

s, 0 , and t are reduced to uncoupled sets of linear 
differential equations in s and t by expanding the variables 
in trigonometric series in the circumferential coordinate 0, 
and treating the nonlinear terms as pseudo loads. The first 
and second derivatives in the meridional coordinate s are 
replaced by the conventional central finite difference 
approximations, ie. 

{z}' = 1/2A ({z} - {z} ) (1) 

i i+1 i-1 



and 




2 

VA ({z}. , 

1+1 



2 Cz} . + £z} . ) 

1 1-1 



(2) 



th 

where {z} is the vector of U, V, W, and M at the i 
i s 



station, A is the uniform dimension between stations, and 

primes denote partial derivatives with respect to s. 
Applying these approximations to the governing set of domain 
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equations leads to 




( 3 ) 



When the shell does not have a pole, fictitious stations 
one increment off of the shell are introduced at each end. 
Both the governing domain equations and the boundary 
conditions are applied at the two boundary points. Thus, 
all finite difference approximations to the derivatives, 
including those of the boundary conditions, are of order 

A 2 

cl . However, prior to the development of SATANS-IxA, the 

treatment of the conditions to be applied at a pole at 
either end of a shell was handled by a simple Euler forward 
or backward difference approximation to the first 
derivative, with truncation error of order A . For example, 
for a pole at s= 0, where i= 1, the first derivative at the 
pole was approximated with 



At the time this procedure for handling the pole 
conditions was developed (1967) it was thought that this 
would not significantly alter the solution. However, it has 
since been discovered that such is not the case. 

For the new pole routine, an expanded forward difference 

^ 2 

approximation of order A is used at s= 0 which takes into 
account the two stations after the pole, instead of just one 
station after the pole as in the Euler scheme. This 
approximation is 



[zf^ = V-A ({2}^ - {Z}^) 



{ 2 } 



1/2^ (-3{Z) + 4{z} 




( 5 ) 



1 



1 



2 
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The conditions to be imposed upon the dependent 
variables at a pole are derived in Reference 14. They are : 

For N= 0, u =v =w‘' =m'' = 0. 

1 1 1 si 



Applying equation (5 ) , these conditions can be put into the 
matrix form 



10 0 0 
0 10 0 
0 0-3 0 
0 0 0-3 




0 0 0 0 
0 0 0 0 
0 0 4 0 
0 0 0 4 




0 0 0 0 

0 0 0 0 

0 0-1 0 

0 0 0-1 




0 . 



where the above 3 matrices are DL, DG, and DF within the 
SATANS programs. 



For N = 



u ±v = d'=w = m = 0, 
1 1 s 



where the plus sign applies at an initial pole, and the 
minus sign at a final pole. The matrix form for these 
conditions is 



-3 0 0 0 
1±1 0 0 
0 0 10 
0 0 0 1 




4 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 




-10 0 0 
0 0 0 0 

0 0 0 0 

0 0 0 0 




For N=2, 

the matrix form is 



u= v= w= m' = 0 
s 



"l 0 

0 1 

0 0 

0 0 




0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 4 




0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0-1 




0 . 



For N > 2, 



u= v= w= 



m = 0 
s 
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and DL= identity matrix, DG= DP= null matrices. 



The solution procedure in 
scheme and starts with 

{z} - [P]^ {z}^ + {X} 

where the values in [P] 

are defined in Reference 

defines a new fP] . 

1 

simultaneously solving the 



SATANS is an elimination 

( 6 ) 

upon the Euler approximation 

The higher order approximation 

This new fP] is obtained by 
1 

pole conditions 



1 

based 
14. 



[DL] {z}^ ♦ [DG] {Z}^ + [DF] = {0}, (7) 

and the domain equation at station 2 next to the pole 

[C] {z} + [B] fz} + [A] £z} = £g} , (8) 

2 1 2 2 2 3 2 

to eliminate {z} Thus, 

£z} = [A]*' (£g} - [C] {z} - [B] {z} ). (9) 

3 2 2 2 1 2 2 

Substituting equation (9) into equation (7) gives 

[DL] [z] + [DG] {z) + [DF] [A]^ ( (g) - C C ] [z] - 

[B] (z) ) = 0. (10) 

2 2 

Combining like coefficients of the [z] vector leads to 

( [DL] - [DF] [Ar[C] ) [Z] + ( [DG] - [DF] [A]'*[B] ) 

2 2 1 2 2 

[z] = - [DF] [A]*' [g] . (11) 

2 2 2 

Finally, solving for [z] yields 
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£Z) = - [DL - DF X A*' X C ] [ DG - DF X A*' X B 1 {z} 

1 2 2 2 2 2 

+ [ DL - DF X A'^ X C^] \ - DF X A*^' ] {g}^. (12) 



-1 

Thus, [P] = - [DL - DF X A*' X C ] [ DG - DF x A'* x B ] and 

1 2 2 2 2 



{X} = [ DL - DF X A"' X C ] [ - DF X A*'] £g} . The new 

1 2 2 2 2 

£ P] matrix has been placed into the "PMATRX" subroutine of 



SATANS-IIA and the new £x) vector has been placed in the 



•'FORCE" subroutine. 



A listing of the pole routine may be found in Appendix 
D. To incorporate this new routine into a SATANS-I or-II 
program, first proceed to the "PMATRX" subroutine and remove 
the fifteen cards that are between, but not including, 
"IF(NN.GT.2) GO TO 90" and "11 CONTINUE".. These cards are 
located after statement number "14" and just before 
statement number "11". Replace the cards removed by the 
ones listed in Appendix D which read from" C IN PMATRX " to 
" 90 M3=MN". Then proceed to the "FORCE" subroutine and 
remove statement number "10". Replace statement number "10" 
with the nine cards listed in Appendix D which read from " C 
IN FORCE" to "DO 11 1= 1,4" . Also place " COMMON 
/IBL5/IBCINL, IBCFNL " into the common area of the " FORCE " 
subroutine. 

This completes the implementation of the new pole 
routine into either SATANS-I or II. 
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IV. PBOB LEM DESCRIPTION 



The geometry of the shallow spherical shell used in this 
study is identical to that used in Reference 1. Briefly, 
the shallow shell can be specified by the non-dimensional 
parameter X , where 

, 2 1/4 1/2 

A = 2[ 3 ( 1 - V ) ] ( H/ h ) . (1) 

H is the rise of the shell, h is the thickness, and V is 
Poisson's ratio. The mass density of the shell is m. All 
shells analyzed had the following dimensions; 



Radii of Curvature 
Thickness 

Modulus of Elasticity 
Poisson's Ratio 



R = R,a = 250 inches 
s ® 

h = 0.25 inches 
E = 30,000,000 psi 
y = 0.3 



All buckling pressures obtained will be listed as a 

percent of the classical buckling pressure of a complete 

sphere, g , where 
0 

2 2 1/2 

q = [ 2 E ( h/R ) ] / [ 3 ( 1 - V ) ] (2) 

0 s 



Forty stations were used over the meridian, 
nondimensional time increment S t, where 



2 1/2 
t= T / ( R^m / E) 
s 



The 



(3) 
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was taken as 0.05 for 3000 time steps, which is a total 
nondimensional time of 150. In addition, the axisymmetric 
analysis was repeated with a larget time step of S t= 0.2 
for a total time of 600. In this study m was selected such 
that t is equal to T. The necessity for the long response 
time is explained in Reference 6. 



In the axisymmetric analysis only the N= 0 harmonic is 
considered. However, in the asymmetric analysis a second 
harmonic is excited by applying an incremental load in that 
harmonic. In addition, analyses of the shells X = 6, 7.5, 
and 11 were made using five harmonics. The step pressure 
load for the axisymmetric harmonic is 



{q 




P (1), 



(«) 



and the step pressure load for the asymmetric second 
harmonic is 



(n) (n) 

{q } = P q^ e {1} , (5) 

(n) 

where n> 0, and f is taken as 0.0001. The value taken 

for the second harmonic in the asymmetric analysis was the 
same as the critical harmonic for the static buckling 
analysis presented by Stilwell and Ball [2]. When there was 
an uncertainty as to which was the critical static harmonic 
the two harmonics in question were both tested. Sun times 
using SATANS-IIA with a two-harmonic analysis for 3000 time 
steps and 40 stations on the meridian took an average of 28 
minutes on the IBM 360/67. 



The parameter used to determine the minimum load at 
which dynamic buckling occurs is the peak value of V, called 



21 



V , where ^ is defined as 
MAX 






( 0 ) 

r H dr / 






dr 



( 6 ) 



r is the normal distance from the axis to the shell, r is 

0 

( 0 ) 

the maximum value of r, M is the normal displacement of 



the axisymmetric response and ^ is the vertical distance 

from the base plane to the undeformed shell. The ? is a 
measure of the volume of the shell deformation. The Fortran 

statements computing V and V are given in Appendix E. 

MAX 

When working a problem that requires these calculations the 

nineteen cards are inserted directly into the "DYNAHIC' 
subroutine right after the *«if” statement that calls the 
"OUTPUT" subroutine. 

For convenience, the response in each asymmetric 

( 0 ) 

harmonic is also measured using equation (6) , with W 

(n) 

replaced with W . The parameter 7 for the asymmetric 



harmonics does not represent a volume of deformation as it 

does for the axisymmetric harmonic. It can, however, be 
used to indicate the relative excitation of the asymmetric 
harmonics. 



The buckling criterion for both the axisymmetric and the 
asymmetric dynamic buckling analysis defines the critical 
load as that load P where a very small increase in P causes 

a very large increase in V . This is the same criterion 

MAX 
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as that used in Ref 



[ 1 ] 
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7. 



RESULTS AND DISCUSSION 



A. STATIC AXISYMHETRIC BUCKLING ANALYSIS 



Table I presents the new results from the static 
axisymmetric buckling analyses for X = 4 through 13 using 
the new pole routine. The two upper curves in Figure 1 
present a comparison of the new results obtained by 
SATANS-IIA with those obtained by Stilwell and Ball [2] 
using the SATANS-I program. As can be seen in this figure, 
fairly significant changes in the buckling load occurred in 
the neighborhood of X = 4,5, and 9; and somewhat smaller 
differences occurred in the region X = 10 through 13. The 
upper data points in Figure 2 present the comparison of the 
new results from SATANS-IIA with those obtained by Huang 
[4]. This comparison shows a very good agreement between 
the two sets of results, except for the largest values of X. 
The new results have eliminated the differences that existed 
between the SATANS-I results and Huang's results. 



B. DYNAMIC AXISYMMETRIC BUCKLING ANALYSIS 

Figure 3 presents the new results for the peak value of 

7 versus P for the various values of X tested. Table 
MAX 

II presents all of the new results for the dynamic 

axisymmetric buckling load. These loads are selected from 
figures constructed just like Figure 3. In every case. 
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except for ^ = 4, a value of P slightly above the P 

CRIT 

value caused a V indicative of buckling, as well as a 

MAX 

nonconvergence of the iterative solution procedure. 



The lower two curves of Figure 1 present a comparison 
versus X of the new axisymmetric dynamic buckling results 
with the previous buckling results obtained by Ball and Burt 
[ 1 ]. In every case the new critical pressure is lower than 
the critical pressure obtained using the Euler approximation 
at the pole. 

The lower data points of Figure 2 present a comparison 
of the new results with those obtained by Huang [5], by 
Stephens and Fulton [6], and by Stricklin £8]. Just as in 
the case of the static axisymmetric buckling analysis, the 
new results compare much more favorably with the other 
results than did the results of Reference 1. It's 
interesting to note that the new results now tend to be 
slightly lower than the other results, whereas the results 
of Reference 1 were higher for almost all values of A . 



C, DYNAMIC ASYMMETRIC BUCKLING ANALYSIS 



Table III presents the new results for the critical 
pressures obtained from the dynamic asymmetric analysis. 
The second harmonics, or critical static harmonics, used in 
the analyses are also presented in Table III. A comparison 
of the critical pressures from the asymmetric analyses. 
Table III, with the critical pressures from the axisymmetric 
analyses. Table II, reveals that only the shell A =6 
buckled at a load below the axisymmetric buckling load. For 
the shell X = 7 the critical buckling load was slightly 
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larger when asymnietric motion was considered. In all other 
cases the buckling was not influenced by the presence of the 
second harmonic. These new buckling results and those by 
Ball and Burt [2] are plotted in Figure 4. The new results 
can be seen to be significantly different from the SATANS-I 
results, where the asymmetric buckling loads were lower than 
the axisymmetric loads for five out of the ten values of 
tested. 

Except for \ = 6 and 7, the relationship between V 

MAX 

and P for the N= 0 harmonic, in the two-harmonic analyses, 

was found to be essentially identical to the relationship 

found in the axisymmetric buckling analysis shown in Figure 

3. Table IV A presents the V versus P data for both the 

MAX 

N= 0 harmonic and the second harmonic, for all values of A 

tested, except for X = 6. Note that, except for X =7 and 

11, V for the asymmetric harmonic is generally very 
MAX 

small, even when the V for the N= 0 harmonic indicates 

MAX 

that the shell has buckled. Thus, except for the shells X = 

6 and 7, the presence of the asymmetric motion does not 
influence the axisymmetric motion, and except for the shells 
X = 6, 7 and 11 the asymmetric motion is very small prior 

to buckling in the axisymmetric harmonic. 

A more detailed analysis of the shell X = 6 has been 

conducted since it was the only shell that revealed any 

significant axisymmetric sensitivity to asymmetric motion. 

This shell was studied using two two-harmonic analyses (N= 

0, 1 and N= 0, 2) and a five-harmonic analysis (N= 0, 1, 2, 

3, and 4) . Figure 5 and Tables IV B and IV C contain values 

of V versus P for both of the asymmetric harmonics, N= 1 
MAX 
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and N= 2 , in the two two-harmonic analyses, as well as the 

values of V for the axisyrametric harmonic, N= 0. Figure 
hAX 

6 and Table IV D present the values of V versus P for the 

HAX 

N= 0, 1,2,3, and 4 harmonics from the five-harmonic study. A 
comparison of the critical buckling load predicted from the 
results of the two two-harmonic analyses in Figure 5 with 
the critical load from the five-harmonic analysis obtained 
from Figure 6 shows that the presence of the additional 
harmonics results in the shell buckling at a slightly lower 
load (0.50), with significant motion in the N= 1 harmonic 
instead of the N= 2 harmonic (see the nocconvergsd solution 
at P= 0.51), which is the critical harmonic for static 
asymmetric buckling. Studies using five harmonics have also 
been conducted for X = 7.5 and X= 11. As can be seen in 
Table IV D the critical harmonic for X = 7.5 remained N= 3; 
however, significant motion occurred in that harmonic at P= 
.41 and .44. In the case of X = 11, relatively large 

asymmetric motion occurred in the asymmetric mode of N= 5 
vice 6 at a value of P= .46. 

The comparison of the new results for the critical 

pressure for dynamic asymmetric buckling with those obtained 

analytically by Stricklin [8], by Akkas [9], and 

experimentally by Lock et al [7] is illustrated in Figure 7. 

The comparison reveals an agreement with Stricklin in every 

case, in general a higher value of P than those obtained 

CHIT 

by Akkas, and most importantly a very good agreement with 
Lock's experimental results. 

When making the comparison between the new results and 
those obtained by Akkas, it is necessary to look at the 
differences in the problem solution parameters used in the 
two studies. For example, buckling results obtained from 
SATANS-IIA using the same time increment as used by Akkas, 
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S t= . 2 for 3000 time steps, were significantly higher 
than those using the time step of 5 t= .05 for many values 
of A . Furthermore, the new results had , in some cases, 
instances of buckling occurring as far out in time as 130. 
Akkas, to shorten computer run times, observed the cap only 
for a time of less than 5. Furthermore, only the harmonics 
N= 1 or 2 or 3 were studied by Akkas for shells ^ = 5 
through 12. If the critical harmonic is not studied, the 
predicted load will be too high. Thus, it appears that 
Akkas' lower bound loads may not be true lower bounds. 

Two additional features of the shell response should be 

noted. First, shells X = 6, 7.5, and 11 exhibited a 

non-buckled response in the axisymmetric harmonic to a load 

larger than the defined critical buckling load. This can be 

seen in Tables IV A and IV C. Second, and most importantly, 

the buckling load proposed by Ball and Burt [1], and used 

here, defines buckling to occur when the V in the 

HAX 

axisymmetric harmonic undergoes a large change due to a 

small change in P. Another criterion for dynamic buckling 

in the asymmetric analysis discussed in Reference 1 is to 

define the buckling load as that threshold load that 

initiates significant growth in the asymmetric harmonic. 

Re-examination of the V versus P data in Table IV A 

MAX 

through D reveals that shells X = 6, 7, and 11 exhibited 

relatively large asymmetric motion at loads smaller that the 

defined buckling load when compared with other V values 

MAX 

for those shells, even though the numbers themselves were 

small when compared with the axisymmetric harmonic. Shells^ 
s7.5 and 12 appear to be borderline cases. If the alternate 
criterion for buckling is used, the critical buckling loads 
for shells X = 6, 7, and 11 become 0.47, 0.45, and 0.45, 

respectively. The shells X = 7.5 and 12 could have buckling 
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loads as low as 0.40 and 0.44, respectively. These values 
are more conservative than the definition based upon 
axisymmetric response. These five shells are the same five 
shells that exhibited an asymmetric buckling load lower than 
the axisymmetric buckling load in Reference 1. 
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VI. 



SUMMASI AND CON CLDSIONS 



A digital computer program for the geometrically 
nonlinear, analysis of totally arbitrarily loaded shells of 
revolution (SATANS-II) was modified to more accurately 
account for the conditions at the pole of the shell. This 
program, called SATANS-IIA, was used to determine the 
buckling load of shallow spherical shells of various sizes 
when subjected to static axisymmetric, dynamic axisymmetric, 
and dynamic nearly axisymmetric step- pressure loads of 
infinite duration. The cap sizes ranged from X = 4 to 13 
including X = 7.5. A comparison was made between the new 
buckling results with the improved pole handling routine and 
the results that did not have the new pole routine. The 
comparison revealed a significant change in buckling 
pressures, due solely to the change from an order 6. finite 
difference approximation of the first derivatives at the 

2 

pole to an approximation of order A . These new critical 

pressures are in very good agreement with the results from 

other studies of the same spherical shells. This good 
agreement with other results, which came about as a result 
of the modification of the pole handling routine, is a 
strong indication that the manner in which the pole 
condition is handled is vital to the accuracy of the 
solutions obtained. 

In the asymmetric analysis, two harmonics were included 
for most of the shells; the axisymmetric harmonic and one 
asymmetric harmonic. Five-harmonic analyses were conducted 
for three of the shells. Two buckling criteria for the 
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asymmetric atalysis were considered. One defined buckling 
as that threshold load that caused a large increase in a 

deformation parameter, V , in the axisymmetric harmonic. 

MAX 

The other, more conservative than the first, defined 

buckling as that threshold load that caused a large increase 

in the V value for the asymmetric harmonic. Both values 
MAX 

have been presented. 



The new static axisymmetric, dynamic axisymmetric, and 
even the dynamic asymmetric critical buckling pressure loads 
appear to be fairly reliable results for perfect, shallow 
shells. The effect of realistic imperfections remains to be 
determined. 



31 



m 



m 

ui 

a: 



O 



Ln 



zn 

in 

LU 



UJ 



O 




H 

I 

CO 

:z 

<c 

H 

<< 

(TJ 



CO 

ra 

CO 

a 

w 



H 

I 

to 

;z 

»< 

H 

< 

CO 



o 

H 

« 

W 

SI 

s 

to 

H 

X . 
-33 



ts 

cn 



J— 

nn 

tx 






SQ 

Ln 






32 




O- 



33 



Figure 2 - CRITICAL STEP-PRESSURE LOAD VERSUS^ 
AXISYMMETRIC (SATANS-IIA VERSUS ALL OTHERS) 
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Figure 3 • PEAK DEFLECTION VERSUS P, AXISYMMETRIC AND 
ASYMaETRIC CASES FOR VARIOUS VALUES OF A (SAIANS-IIA) 
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Figure 5 - PEAK DEFLECTION VERSUS P FOR THE ASYMMETRIC 
ANALYSES OF A = 6 (N=0,1 AND N=0,2) 
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Figure 7 - CRITICAL STEP-PRESSURE LOAD VERSUS A 
ASYMMETRIC ANALYSES (SATANS-IIA VERSUS ALL OTHERS) 



A 



TABLES 



1 . TABLE I Critical press u re loa ds from the static 
Misisaetric an alyses . 



A 


4 


5 


6 


7 


8 


9 


10 


1 1 


12 


13 


P 

CRIT 


.568 


.616 


1.0 


1.048 


1.12 


.936 


. 83 2 


. 832 


1.016 


1.032 



2 . TABLE II Critic al step-press ure lo ad s £r om the 
ax i symmetric d ynamic analy se s . 



A 


4 


5 


6 


7 


7.5 


8 


9 


10 


1 1 


12 


13 


P 

CRIT 


.45 


.44 


.59 


.53 


. 42 


. 40 


.39 


.33 


.47 


.45 


. 35 



3. TABLE ^I C ritical step-press ure loa ds from the 
dy n ami c a symm e tri c analyses and critical asy m metric 
harmonics. 



A 


5 


6 


7 


7.5 


8 


9 


10 


1 1 


12 


13 


P 

CRIT 


.44 


.52 


.54 


. 42 


.40 


.39 


.33 


.47 


.45 


.35 


N 

CRIT 


1 


2 


3 


3 


4 


5 


6 


7 


8 


9 



39 




I 



I 



I 



4 



versus 



P. 



TAB LE IV Dy namic asymm^ric analy ses for 



1. TABLE IV A. Two-harmonic analyses for all values of 
A exept A = 6 . 

A = 5 N= 0 and 2 N= 0 and 1 



p 


.43 


. 44 


.45 




P 


.44 


. 45 


o 

II 

VZi 


. 1659 


. 1676 


.6606 




N= 0 


. 1675 


.6606 


N= 2 


,.0004787 


.0000566 


.0687 




N= 1 


.0003145 


. 0687 












P 


.46 












o 

II 


.7653 












N= 1 


.001092 



\ = 7, N= 0 and 3 



P 


.45 


.46 


.47 


• 

CD 


.49 


.50 


. 52 


o 

II 


.09452 


.09571 


.09812 


.1005 


. 1029 


.1052 


. 1099 


N= 3 


.000889 


.007456 


.05052 


.04323 


.0279 


.0335 


. 07488 


P 


.53 


. 54 


.55 






N= 0 


. 1122 


.1146 


.2709 




N= 3 


. 05997 


.06252 


.03809 





A = 7.5, N= 0 and 3 



P 


.40 


.4 1 


.42 


.43 


.44 


.45 


N= 0 


.0703 


.07228 


.07429 


.2636 


.07837 


.2076 


N= 3 


.0001094 


.001 296 


.0004304 


.002754 


.001188 


.000338 


P 


. 46 




o 

II 


. 200 


N= 3 


.0003276 



40 



A = 8, N= 0 and 4 



p 


• 

U) 

CD 


.39 


.40 


.41 


.42 


.43 


N= 0 


.05893 


.0607 


.0624 


.1964 


. 1713 


. 1957 


tl 

:z5 


.0000566 


.0000703 


.0000364 


.0000333 


. 0000274 


.0000299 


p 


. 44 




o 

tt 


. 2297 


N= 4 


.0000326 



A = 9 N= 0 and 4 N= 0 and 5 



P 


.38 


.39 


.40 


P 


.40 


o 

It 


.04738 


.04875 


. 1576 


N= 0 


.05012 


II 


.00003597 


.00004635 


. 00004497 


N= 5 


.00008385 



X =10, N=0and5 



P 


.32 


.33 


.34 


.36 


. 38 


.40 


' 2 
II 

O 


.03239 


. 03347 


.1086 


.1217 


. 1288 


. 1235 


N= 5 


.0000281 


.0000472 


.00004125 


.00002103 


. 0000449 


.000114 



A = 11, N= 0 and 6 



P 


. 45 


. 46 


.46 


• 

CD 


.49 


. 50 


N= 0 


. 03910 


.04004 


.04099 


.09814 


. 0424 1 


.08824 


N= 6 


.004595 


.01332 


.02232 


.02864 


.03955 


.02813 



A = 12, N= 0 and 7 



P 


.44 


.45 


.46 


N= 0 


.03236 


.03316 


. 08633 


N= 7 


.00004214 


.0004561 


.00005158 



A = 13, N= 0 and 8 



P 


.34 


.35 


.36 


.38 


. 40 


o 

II 


.02119 


. 02185 


.06637 


.07844 


.07381 


00 

II 


.00001 148 


.00001134 


.000006607 


.000008245 


.000119 
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TABLE IV B. Two-harmonic analyses with N= 0 and 1 , 
A = 6 only. 



p 


.48 


.49 


.50 


.51 


.52 


.53 


N= 0 


. 1350 


. 1385 


.1421 


. 1460 


. 1499 


.4453 


N= 1 


. 0002797 


.000195 


.000245 


.000926 


.04081 


.3668 



3. TABLE IV C. Two-harmonic analyses with N= 0 and 2, 
A = 6 only. 



P 


.44 


.45 


.46 


.47 


.48 


.49 


.50 


o 

II 


.1218 


. 1250 


.1276 


. 1320 


.1350 


. 1385 


. 1479 


CN 

It 


. 000239 


.0101 


.0293 


.01976 


.08768 


. 1223 


.1419 


P 


.51 


.52 


. 53 


.54 


.55 


.56 


o 

II 


. 1526 


.1499 


.5137 


.5060 


.2040 


.5305 


N= 2 


. 1 654 


.1816 


.3878 


.3996 


.2156 


.3617 



4. TABLE IV D Five-harmonic analyses for selected 
shells. 



A = 6 N= 0, 1 ,2,3, and 4 



P 


.47 


• 

00 


.49 


.50 


.51 


.52 


.53 


N= 0 


. 13 13 


. 1347 


.1382 


.1460 


.3797 


.1498 


. 2200 


N= 1 


.00021 


.02108 


.02215 


. 003676 


.3743 


.01276 


.0616 


N= 2 


.0187 


.0953 


.1069 


. 1507 


.1502 


.1279 


. 2385 


ss ' 
II 

UJ 


.000181 


.006237 


.01437 


.00163 


. 04 05 


.0123 


. 03978 


N= 4 


.0031 


.04757 


.05428 


. 04402 


. 05896 


.0495 


. 064 


P 


.54 




N= 0 


.3854 


N= 1 


.3953 


N= 2 


.1671 


N= 3 


.05298 


II 


.0613 
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A = 7.5 N= 0,1, 2, 3,4 and 4 



p 


.40 


.41 


.42 


.43 


.44 


.45 


N= 0 


.0703 


.07228 


.07429 


.2592 


.07837 


.2544 


N= 1 


.00004855 


.00004198 


.00006093 


.0002737 


.0001167 


.005952 


N= 2 


.0001164 


.00007456 


.0004184 


.0000982 


.000788 


.0003188 


N= 3 


.0001277 


.001187 


.0004597 


.0002853 


.00107 


.0003188 


tt 


.0008224 


.0001898 


.0002448 


.0000526 


.000280 


.000134 



A = 11 N= 0,4, 5, 6, and 7 



P 


.45 


.46 


.47 


.48 


.49 


.50 


II 

o 


.03910 


.04004 


.0499 


.04195 


.04291 


.1040 


N= 4 


.0005759 


.001 263 


.0009774 


.001388 


.002657 


.00 1565 


N= 5 


.009568 


.0140 


.0124 


. 02239 


.02759 


.01 548 


N= 6 


.002560 


.007828 


.02330 


. 02767 


. 02602 


.02644 


N= 7 


.0001743 


.0001486 


.0002021 


.01202 


.02048 


.02064 
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APPENDIX A 



LISTING OF SATANS-IIA 
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3 XIQMIHA 



Azi'EUTiiz HOt aaiPd to^ii 



INPUT DATA GUIDE FOR SATANS-It SATANS-II, AND 

SATANS-IIA 



CARD 


CCLLMNS 


FORMAT 


ITEM 


EXAMPLE 


MEANING 




1 


1-72 


18A4 


TITLE 


- 


ENTER ANY 72 


CHARACTERS 


2 


1-5 


15 


NO 


1 


THE PROBLEM 
0<NC<10000. 


NLMBER, 


2 


e-ic 


L5 


SDYNAMC 


F 


FOR A STATIC 


ANALYSIS, 



SET SDYNAMC = F, 

T FOR A DYNANIC ANALYSIS, 

SET $DYNAMC = T. 



2 


11-15 


15 


IMODE 


0 


FCR NO MODAL OUTPUT DATA 










1 


FCR MCDAL OUTPUT DATA. 


2 


16-20 


15 


NDIMEN 


0 


DIMENSIONAL CLTFUT DATA. 










1 


NCNCIMENSICNAL OUTPUT. 


2 


21-25 


15 


NTHMAX 


8 


SUMMED SOLUTION WILL EE 
PRINTED AT NTHMAX MERID- 
IANS, 0<=NTHMAX<=36. 


2 


26-30 


15 


IFREQ 


2 


SOLUTION WILL BE PRINTED 
AT THE FIRST STATION, 
EVERY SUBSEQUENT IFREQ 
STATION ANC THE LAST 
STATION, 0<IFREC<=KMAX. 


2 


21-35 


I 5 


IPRINT 


3 


EVERY IPRINT CONVERGED 
SOLUTION WILL BE PRINT- 
ED. 


2 


36-40 


15 


IBC INL 


-1 


IF THE SHELL HAS A POLE 
AT THE FIRST STATION. 










0 


IF THE SHELL HAS NO POLE 
AT THE FIRST STATION. 


2 


Al-45 


15 


IBCFNL 


-1 


IF THE SHELL HAS A POLE 
AT THE LAST STATION. 










0 


IF THE SHELL HAS NO POLE 
AT THE LAST STATION. 
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CARC CGLLf'NS FORMAT ITEM 
2 46-50 15 KMAX 

2 51-55 15 MNMAX 

2 56-6C 15 MAXM 

2 61-65 15 LSMAX 

2 66-70 15 LCHMAX 

2 71-75 15 ITRMAX 

2 76^80 15 IC 



RLE MEANING 

NUMBER OF MERICICNAL STA- 
TICNS. 

NOTE: KMAX<201 FCR SATANS 
-I WITHOUT PLOTS AND 
KMAX<101 FOR SATANS-I 
WITH PLOTS CR FOR SATANS 
-II. SATANS-IIA IS 
UNLIMITED. 



NUMBER OF SERIES COEFFI- 
CIENTS USED TO DESCRIBE 
THE INITIAL CONDITIONS, 
PRESSURE ANC THERMAL 
LOADS (AND INITIAL IMPER- 
FECTICNS IF USING SATANS 
-II OR IIA). MNMAX<=MAXM. 



MAX NUMBER OF HARMONICS 
IN THE SOLUTION, LIMITED 
TO 99. 



FOR A LINEAR ANALYSIS. 
USE MANY LOAD STEPS FCR 
A NONLINEAR STATIC 
ANALYSIS. 

FOR A DYNAMIC ANALYSIS, 
LSMAX IS THE NLM8ER OF 
TIME INCREMENTS, WHERE 
LSMAX= T MAX/^T. 



THE NUMBER CF LOAC STEP 
SIZE REDUCTIONS IN A 
STATIC ANALYSIS, RECCM- 
MENDED RANGE = 2-4. 

FOR A DYNAMIC ANALYSIS. 



FOR A LINEAR ANALYSIS. 

THE NUMBER CF ITERATIONS 
AT A LOAD OR TIME STEF. 
FOR A NONLINEAR ANALYSIS, 
SUGGESTED RANGE = 10-30, 
UP TO 50 FOR SPECIAL 
CASES. 



INITIAL CONDITIONS- SET 
TO 0 FOR A STATIC ANALY- 
SIS, OR FOR A DYNAMIC 
ANALYSIS WHERE THE SHELL 
IS AT REST AT T= 0. 

FOR A DYNAMIC ANALYSIS 
WITH INITIAL CCNDITICNS. 



EXAM 

35 

7 

7 

1 

99 

3000 

2 

0 

1 

30 

0 

1 
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C^RC 


COLUMNS 


FORMAT 


ITEM 


EXAMPLE 


MEANING 


3 


1-12 


E12.3 


NU 


0.3 


POISSON'S RATIO, V . 


3 


12-24 


512.3 


SIGO 


1000. 0 


REFERENCE STRESS LEVEL, 

. 










1.0 


IF TFE INPUT D^TA IS 
DIMENSIONAL. 


3 


24-36 


E12.3 


ELAST 


.3E8 


REFERENCE MCCULUS OF 
ELASTICITY, E- . 

IF TFE INPUT D^TA IS 










1.0 












DIMENSIONAL. 


3 


37-48 


E12.3 


TKN 


.4E-2 


REFERENCE THICKNESS, . 










1.0 


IF THE INPUT D4TA IS 
DIMENSIONAL. 


3 


4S-6C 


E12.3 


CHAR 


8.16 


CHARACTERISTIC SHELL 
DIMENSION, a . 










1.0 


IF THE INPUT DATA IS 
DIMENSIONAL. 


3 


61-72 


E12.3 


TEEO 


0.0 


IF A STATIC ANALYSIS. 








• 


996E-5 


REFERENCE TIME, . 


4 


1-12 


E12.3 


DELCAD 


0.2 


FOR A STATIC ANALYSIS, 
DELCAD IS THE LCAC INCRE- 
MENT. IT REMAINS UN- 
CHANGED UNTIL THE SOLU- 
TION FAILS TO CCNVERGE IN 
ITERMAX ITERATIONS, UhEN 
IT IS REDUCED BY A FACTOR 
OF FIVE. A MAXIMUM CF 
LCHMAX SUCH REDUCTIONS 
WILL CCCUR. 








• 


L823E-6 


FOR A DYNAMIC ANALYSIS, 
DELCAD IS THE NCNDIMEN- 












SICNAL TIME INCREMENT. 


4 


13-24 


E12.3 


EPS 


0.01 


THE CONVERGENCE CRITERION 
RECOMMENCED RANGE OF 
0.0KEPS<0.CC1 . 


CARC 


4A IS ONLY INCLUDED FOR A SATANS-II OR SATANS-IIA RUN. 


4^ 


1-5 


15 


JUMP 


1 


FOR AN ANALYSIS USING 
SINGLE SERIES EXPANSIONS. 










2 


FOR AN ANALYSIS USING 
DOUBLE SERIES EXPANSIONS. 


4A 


5-lC 


15 


MPERFS 


0 


AN ANALYSIS WITHOUT IM- 








1 


PERFECTIONS. 

AN ANALYSIS WITH IMPERFEC- 





TIONS. 

NOTE: IF JUMF=28 MPERFS 
MAY BE 0 OR 1. IF JU^P 
=lt MPERFS MUST BE 0. IF 
MPERFS=1, JUMP MUST BE 2. 
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CARC CCLtNN FCRMAT ITEf< EXAMPLE MEANING 



INCUDE AS MANY CARDS 5 AS NECESSARY TC SPECIFY NThMAX 
MEPICIANS. IF NTHMAX EQUALS 0, OMIT CARC 5, 

5 1-72 6E12.3 10,0 A LIST OF CIRCUMFERENTIAL 

COORDINATES d ♦ IN DEGREES 
ANC TENTHSt WHERE THE 
SCLUTICN PRINTCUT IS CE- 
SIREC. THE LIST MUST HAVE 
NTHMAX ENTRIES. 



IF IBCINL= -It CMIT CAROS 6 THROUGH 14. IF IBCFNL= -1, CMIT 
CARDS 15 THROUGH 23. CARDS 6 THROUGH 23 DESCRIBE THE 
BCUNCARV CCNDITICNS AT THE FIRST, ANC THEN AT THE LAST STA- 
TION. THE BOUNDARY CONDITIONS EXIST CN THE TOTAL VARIABLES, 
NOT CN THE INDIVIDUAL HARMONICS. LOADINGS APPLIED THROUGH 
SPECIFICATION CF BOUNDARY CCNOITIONS ARE TAKEN IN THE ZEPC- 
ETH HARMCNIC (N=C) ONLY, AS THE COLUMN MATRIX f4] IS SET TC 
ZERC FOP HARMONICS GREATER THAN ZERO. THE BOUNDARY CON- 
OITICNS ARE DIMENSIONAL. THE FORMAT CF CAROS 6 TFRCUGH 23 
IS AE16.8. 



CARC 6rl5 


CARD 7,16 


CARD 8,17 


CARD S, 


18 




A ( 1 , 1 ) 


A (1,2) 


A (1,3) 


A (1,4)“ 




N, “ 


A (2,1) 


A (2,2) 


A (2,3) 


A (2,4) 




N/« 


A ( 3 ,1) 


A (3,2) 


A (3,3) 


A (3,4) 






A(^,l) 


A (4,2) 


A (A, 3) 


A (4,4) 




U ^ J 



CARC 10 
HA ( 1,1 ) 
yt(2,l) 
a(3,1) 
>U4,1) 



IS CARD 11,20 
A (1,2) 
( 2 , 2 ) 
A.(3,2) 

A (4,2) 



CARD 12,21 
A (1,3) 

A (2,3) 

A (3.3) 

A (4,3) 



CARC 13,22 ^ 
yV(l,4) U 

A (2,4) V 

A (3,4) W 

A ( 4 , 4 ) M , 



CARD 14 ., 23 
^ (ir 
= T (2) 

Ji (3 ) 

/(4) 



CARC 24 IS: 

1. INCLUDED FOR A SATANS-I STATIC ANALYSIS. 

2. INCLUDED BUT BLANK FCR A SATANS-I DYNAMIC ANALYSIS. 

3. CMITTEC FCR A SATANS-II ANALYSIS. 

A. INCLUDED BLANK FOR DYNAMIC USED FCR STATIC SATANS-IIA 
ANALYSES. 



CARC COLUMN FORMAT ITEM EXAMPLE 
24 1-2 L2 $PLOTS F 

T 



24 3-4 



L2 $MCDAL F 

T 



MEANING 

INDICATES PLOTS ARE NCT 
DESIRED. 

INDICATES PLOTS ARE 
DESIRED; 

INCICATES PLOTS ARE FOP 
SUMMED SOLUTIONS ONLY. 
INCICATES PLOTS ARE FCR 
MODAL SOLUTIONS ONLY. 
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FOR TFE REMAINDER GF CARO 24 ENTRIES, C INDICATES THAT NO 
PLOTS APE DESIRED FOR THE PARTICULAR ITEM, AND 1 INDICATES 
THAI THEY ARE DESIRED. ALL GRAPHS ARE PLOTTED AS THE IN- 
DICATED ITEM VERSUS THE STATION NUMBER. IF A COMPLETE PLCT 
IS DESIRED, INSUTE IFREC = 1. 



CARD 


COLUMN 


FORMAT 


ITEM EXAMPLE 


MEANING 


24 


5- 6 


12 


IRAOI I 


1 


PLCT THE RADII AS COMPUT- 
ED BY SUBROUTINE GEOM. 


24 


7- 8 


12 


IGAMMA 


1 


PLCT AS CCMFUTED EY 

SUBROUTINE GEOM. 


24 


S-IC 


12 


lOHEGS 


1 


PLCT AS COMPUTED EY 

SUBROUTINE GEOM. 


24 


11-12 


12 


lOMEGT 


1 


PLCT AS COMPUTED EY 

SUBROUTINE GEOM. 


24 


12-14 


12 


IDEOMS 


1 


PLCT (*>'$ AS COMPUTED EY 
SUBROUTINE GEOM. 


24 


15-U 


12 


IBSTIF 


1 


PLCT THE. STIFFNESS D AS 
COMPUTED BY SLERCUTINE 



BDE. 



24 


17-18 


12 


IDSTIF 


1 


PLOT THE STIFFNESS D AS 
CCMFUTED BY THE SUBROU- 
TINE BOB. 


24 


lS-20 


12 


IBBSTf 


1 


PLCT THE STIFFNESS 

AS COMPUTED EY SLERCUTINE 

BCE. 


24 


21-22 


12 


IDDSTF 


1 


PLCT THE STIFFNESS dd/dS 
AS COMPUTED EY SUBROU- 
TINE BOB. 


24 


23-24 


12 


IPR 


1 


PLOT THE NORMAL COMPONENT 
OF THE PRESSURE LOAD. 


24 


25-26 


12 


IPS 


1 


PLOT THE MERIDIONAL CCM- 
PCNENT OF THE PRESSURE 
LOAD. 


24 


27-28 


12 


IPT 


1 


PLCT THE CIRCUMFERENTIAL 
COMPONENT OF THE PRESSURE 
LOAD. 


24 


2S-2C 


12 


ITT 


1 


PLCT THE THERMAL LCAD. 


24 


31-22 


12 


IMT 


1 


PLOT THE THERMAL MOMENT. 


24 


33-24 


12 


IDTT 


1 


PLCT d/dSoF the THERMAL 
LOAD. 


24 


25-26 


12 


IDMT 


1 


PLOT d/cHf OF THE THERMAL 
MOMENT. 


24 


27-38 


12 


INS 


1 


PLOT THE MERIDIONAL MEM- 
BRANE FORCE DISTRIBUTICN. 
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CARC 


COLUMN FORMAT 


ITEM EXAMPLE 


MEANING 


24 


39-4C 


12 


INTH 


1 


PLOT THE CIRCUMFERENTIAL 
MEMBRANE FORCE DISTRIBU- 
TION. 


24 


^1-42 


12 


INSTH 


1 


PLOT THE MERIDIQ- 
CIRCUMFERENTI AL MEMBRANE 
FORCE DISTRIBUTION. 


24 


43-44 


12 


IQS 


1 


PLOT THE TRANSVERSE FORCE 
DISTRIBUTION. 


24 


45-46 


12 


IMS 


1 


PLOT THE MERICICNAL MOM- 
ENT DISTRIBUTION. 


24 


41-4E 


12 


IMTH 


1 


PLOT THE CIRCUMFERENTIAL 
MOMENT DISTRIBUTION. 


24 


49-50 


12 


IMSTH 


1 


PLOT THE MERICIC- 
CIRCUMFERENTIAL MOMENT 
DISTR IBUTION. 


24 


51-52 


12 


lU 


1 


PLOT THE MERIDIONAL DIS- 
PLACEMENT DISTRIBUTION. 


24 


53-54 


12 


IV 


1 


PLOT THE CIRCUMFERENTIAL 
DISPLACEMENT DISTRI- 
BUTION. 


24 


55-56 


12 


IW 


1 


PLOT THE NORMAL DISPLACE- 
MENT DISTRIBUTION. 


24 


57-56 


12 


IPHIS 


1 


PLOT THE MERICICNAL ROTA- 
TION DISTRIBUTION. 


24 


59-60 


12 


IPHIT 


1 


PLOT THE CIRCUMFERENTIAL 
ROTATION DISTRIBUTION. 


24 


61-62 


12 


IPHI 


1 


PLOT THE MERIDIC- 
CIRCUMFERENTIAL ROTATION 
DISTRIBUTION. 


INCERT IMPERFECTION 


DATA HERE 


FOR 


A SATANS-II OR SATANS-IIA 


ANALYSIS ITH 


IMPERFECT IONS. 


INSURE FORMAT OF THE IMPERFEC- 


TICN 


DATA IS 


COMPATIBLE WITH 


THAT 


SPECIFIED IN THE USER- 


WRITTEN SLBRCLiTINE 


IMPERF. 






25 


1-2 


12 


IRNAGN 


0 


INDICATES THES IS THE 
ONLY RUN. 










1 


INDICATES ANOTHER RUN IS 



TO BE MACE. ACC ANOThEP 
COMPLETE SET CF DATA 
CARCS AFTER THIS CARD IS 
IRNAGN= 1. 
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APPENDIX D 



LISTING OF NEW POLE ROOTINE FOR SATANS-IIA 
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APPENDIX E 



LISTING 



OF CARDS FOR V AND V 

MAX 
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ctfA ? iiO\ 20S40 to Mmur 



THE FCLLCVvING CAROS ARE TC BE PLACED I^TC THE CYKAMIC SUBRCLTINE IF NEECEO 



>T 

• 

O 

rM 

UJ 

\T\ 



CJ 

o 



Qd 

o 



Qd 

<1 

CD 

> 

to 



Qd 




z- 


CD 


< 




X •• 


< 


OJ 




<x>r 


II 


> 








>T 




•k •• 


X 








w 






II 


OQ 






XQ. 


> 


-JX 




**QJ 


< 


Z)<t 

os: 


CO 


XOO 






o 


w 




<w 


O' 


QCUJ 


• X 


Odd 




<x 


Ow 




^ a 


CD*-* 


II CD 




^ 1- 


>»- 




H-X 


w 


w 


X< 


< 


Qd O 




w • 


zs: 


-Rr S: O 




CDI- 




— UJUJ 


UJ 


>o 




h-OZ -- 


h-Qd 


< • 


x*-«z 


Xt/^UJ^O 


00< 


X 


XX 


eg -«• Q-f ^ 


OCQ 


X— *— 


Qds:< 


XO CO-J^H- • 


W> 


<-iX 


oos: 


wujicyiiu 


— ^ m 


X 


U-* - 


2. MOOUJ^ 


rr) •* 


•‘CJQd 




it:rg + zi- • 


COX 


^UJ<X 


OOCM tl 


•R- II xacxof-i— 


rHin 


It tcQ 


1—^2. 




Z • 


• -« +ZZ— II UJ II >Ow^ 


UJw 



UU tl >oo 

2:^00 II 

1 —^ -I 
<IUJLJ^ 
H-LJUICJ 
CO 

U\ 

cu 



2100+ II II ^«/>h-<>OUJ<l»-»a0(/>0D 
W ,H 2. a. QC UJ UJ 2. ^ ^ <t 

II II ww 

^ W-jr— ^C.UJ»~U-r-tiLJ LJC-M-L.U* 

^ LJ it Ul CJ -> h-ii-h 1-^ j; UL o UJ H-I •— * 



sr 

cu 



>u 



0)01 P*- 

C0^4J GO 



l62 



LIST OP FIGURES AND TABLES 



A. FIGURES 

1. Critical step-pressure load versus 1 axisymmetric 

analyses (SATANS-I versus SATANS-IIA) 32 

2. Critical step-pressure load versus X axisymmetric 

analyses (SATANS-IIA versus all others) 33 

3. Peak deflection versus P, axisymmetric and asymmetric 

cases for various values of X 34 

4. Critical step-pressure load versus X asymmetric 

analyses (SATANS-I versus SATANS-IIA) 35 

5. Peak deflection versus P for the asymmetric analyses of 

X = 6 (N=0,1 and N=0,2) 36 

6. Peak deflection versus P for the asymmetric analyses of 

X = 6 (N=0, 1,2,3/and4, only N=0,1,and2 plotted) 37 

7. Critical step-pressure load versus A asymmetric 

analyses (SATANS-IIA versus all others) 38 



B. TABLES 

1. Critical pressure loads from the static axisymmetric 

analyses 39 

2. Critical step-pressure loads from the axisymmetric 



163 



/ 



dynamic analyses ^ 39 

3. Critical step-pressure loads from the dynamic asym- 
metric analyses and critical asymmetric harmonics... 39 

4. Dynamic asymmetric analyses for V versus P 40 

MAX 



164 



LIST OF REFERENCES 



1. Ball, R. E. and Burt, J. A., " Dynamic buckling of 

shallow spherical shells". J ournal of A pplied 
Mech anics , v. Paper No. 73-APM-A, p. 411 to 416, 
June 1973. 

2. Stilwell, W. C., and Ball, R. E., "A Digital Computer 
Study of the Buckling of Shallow Spherical Caps and 

Truncated Hemispheres, " NASA CR-1998, June 1972. 

3. Ball R. E., "A Program for the Nonlinear Static and 

Dynamic Analysis of Arbitrarily Loaded Shells of 

Revolution," Cump u ters and Str uctures, Vol. 2,1972, p. 
141 to 162; also "A Computer Program for the 

Geometrically Nonlinear Static and Dynamic Analysis of 
Arbitrarily Loaded Shells of Revolution, Theory, and 
Users Manual," NASA CR-1987, April 1972. 

4. Huang, N. C., " Unsymmetrical Buckling of Thin Shallow 
Spherical Shells", J our nal of A pplied Mechanics , v. 
31, p. 447 to 457, September 1964. 

5. Huang, N. C., " Azisymmetric Dynamic Snap-Through of 

Elastic Clamped Shallow Spherical Caps Due to Centrally 
Distributed Pressures", MAA Journal , v. 7, No. 2, p. 
215 to 220, Feb. 1969. 

6. Stephens, H. B., and Fulton, R. E. , " Azisymmetric 
Static and Dynamic Buckling of Spherical Caps Due to 
Centrally Distributed Pressures", AIAA Journ al, v. 7, 
No. 11, p. 2120 to 2126, November 1969. 

7. Lock, M. H. , Okubo, S. and Whittier, J. S. , 



165 



p: 




"Experiments on the Snapping of a Shallow Dome under a 
Step Pressure Load", A IAA J ourn al , v. 6, No. 1, p. 

1320 to 1326, July 1968. 

8. Stricklin, J. A., et al. , " Nonlinear Dynamic Analysis 
of Shells of Revolution by Matrix Displacement Method", 
AIAA Journal, v. 9, No. 4, p. 629 to 636, April 
1971. 

9. Akkas, N., ” Bifercation and Snap-Through Phenomena in 

Asymmetric Dynamic Analysis of Shallow Spherical 
Shells", Co mp u te rs and Structure s, v. 6, p. 241 to 

251, March 1976. 

10. Ryan, E. A., "A Digital Computer Study of the Buckling 
of Actual Imperfect Cylinders - A Modification of the 
Computer Program SATANS - Theory and User's Manual for 
SATANS- I and SATANS-II," A. S. Thesis, Naval 
Postgraduate School, June 1970. 

11. Ball, H. E. and Ryan, B. A., " Computer Analysis of 

Buckling of Imperfect Shells", Journa l of the 

Str uctural Division , v. 99, No. ST10, p. 2097 to 
2108, October 1973. 

12. Ball, R. E. et al. , " A Comparison of Computer Results 
for the Dynamic Response of the LMSC Truncated Cone", 
C om pute rs and Structure s, v. 4, p. 485 to 498, July 
1973. 

13. Aerospace Structures Information and Analysis Center, 
APFDL/FBR, 

Wright-Patterson AFB, 

Ohio 45433 



166 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 

3. Department Chairman, Code 67 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 

4. Assoc Professor R. E. Ball, Code 67Bp 
Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

5. Professor Gilles Cantin, Code 69Ci 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

6. Professor Johann Arbocz 
Firestone Flight Science Laboratory 
California Institute of Technology 
Pasadena, CA. 91109 

7. LT Michael D. Shutt, USN 
7186 Emerald Ave. 

Dublin, California 94566 



167 



Copies 

2 

2 

1 

1 

1 

1 



168183 




Thes i s 
S49325 
c .1 



Thesis 

S49325 Shutt 
c.l Static and dynamic 

buckling of shallow 
spherical shells sub- 
jected to axi symmetric 
and nearly axlsymr 
metric step pressure 
loads using SATANS-IIA 
a modified version of 
SATANS-IL. 



163188 

Shutt 

Static and dynamic 
bi'c''1ing <^f shallow 
spherical shells sub- 
jected to axi symmetric 
and nearly axisym- 
metric step pressure 
loads using SATANS-IIA, 
a modified vers -on of 

SATANS-II. 



thesS49325 

dynamic buckling of shallow s 




3 2768 001 95377 1 
DUDLEY KNOX LIBRARY 



